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In  contrast  to  the  rapid  decrease  predicted  by  conventional 

C  ’  theory,  the  magnetostriction  ^OO  ^-ron  has  a  large  maximum  just 

r' 

below  the  Curie  temperature.  We  propose  a  mechanism  based  on  the 
fact  that  an  ellipticity  in  the  quasi-particle  spectrum  permits  a 
lowering  of  the  free  energy  by  distortion;  an  equivalent  mechanism 
|  arises  from  the  anisotropic  magnon-phonon  interaction  near  the  zone 
boundary.  This  latter  interaction  is  large  at  temperatures  such 
]  that  magnon  renormalization  (due  to  magnon-raagnon  interaction) 

!  lowers  the  magnon  spectrum  to  degeneracy  with  phonons  at  the  zone 
edge.  The  degeneracy  temperature  agrees  well  with  the  temperature 

! 

i  of  the  X100  maximum  in  iron.  Adding  silicone  raises  impurity  states 
from  the  phonon  spectrum  and  thence  lowers  the  degeneracy  tempera¬ 
ture,  but  increases  the  range  of  temperature  over  which  near¬ 
degeneracy  occurs;  this  agrees  with  the  observed  shift  and  broaden¬ 
ing  of  the  PeaK*  The  degeneracy  does  not  occur  in  nickel, 

nor  does  the  ^-j.00  Pealc*  The  mechanism  also  predicts  a  monotori cally 
decreasing  of  the  opposite  sign  to  l^oo’  aS  °hserved  in  iron. 


*Supported  in  part  by  the  Office  of  Naval  Research^ 


I.  INTRODUCTION 


The  conventional  magnetoelastic  coupling  theory  of  mag- 

12  3 

netostriction  *  and  the  observations  of  Tatsuraoto  and  Okamoto'' f 

on  iron  are  in  puzzling  disagreement.  The  magnetoelastic  coupl¬ 
ing  theory  predicts  that  the  magnetostriction  of  ferromagnets 
should  fall  monotonically  to  zero  with  increasing  temperature.'. 

In  contrast  the  magnetostriction  constant  ^q0  iron,  shown  $ 
in  Fig.  la,  increases  with  increasing  temperature,  exhibiting  a 
minor  maximum  at  about  one-fifth  of  the  Curie  temperature  and  a 
major,  large  maximum  just  below  the  Curie  temperature.  Addition 
of  small  fractions  of  silicon  shifts  the  major  maximum  to  lower 
temperature  and  broadens  it.  The  second  magnetostriction  con¬ 
stant  shown  in  Figure  lb,  falls  monotonically  to  zero., 

Finally,  the  peak  in  ^00  is  absent  in  nickel. 

In  this  paper  we  propose  a  mechanism  which  accounts 
qualitatively  for  the  magnitude  and  location  of  the  high-tempera- 
ture  maximum  of  in  * ron *  for  shift  and  broadening  of  the 

maximum  with  the  addition  of  silicon,  for  the  absence  of  the  ef¬ 
fect  in  nickel,  and  provides  criteria  for  the  presence  of  the 
effect  in  other  materials.  Furthermore,  the  theory  properly  pre¬ 
dicts  the  monotonic  behavior  of  and  the  fact  that  an<* 

^100  are  opposite  sign. 

The  mechanism  depends  on  the  existence  of  an  asymmetry  in* 
the  excitation  spectrum  of  the  system,  this  asymmetry  being  deter¬ 
mined  by  the  axis  of  the  magnetization.  For  clarity  we  first 
illustrate  the  effect  by  assuming  a  simple  eccentricity  of  the 
spin-excitation  (raagnon)  spectrum,  such  as  results,  for  instance, 

k 

from  the  dipolar  interaction  of  the  spins. 


f 
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Such  an  eccentricity  is  indicated  schematically  in  Fig.  2, 
where  the  surfaces  of  constant  frequency  of  the  excitations  are 
shown  as  elliptical. 

The  free  energy  of  the  crystal  is  the  sum  of  individual 

contributions  FGc^T)  from  the  individual  modes.  The  product  of 

F(l?,T)  and  of  the  density  of  modes  in  reciprocal  space,  *rs  , 

(2n)-J 

defines  a  scalar  field  in  reciprocal  space.  This  free  energy 
density  must  be  integrated  over  the  Brillouin  zone  to  compute 
the  free  energy  of  the  crystal. 

The  surfaces  of  constant  free  energy  density  coincide 
with  the  surfaces  of  constant  frequency  of  the  modes,  and  an 
eccentricity  of  the  latter  implies  a  similar  eccentricity  of  the 
former. 

If  the  crystal  is  distorted,  two  effects  must  be ' considered. 
First,  the  boundaries  of  the  Brillouin  . zone  over  which  the  free 
energy  density  is  to  be  integrated,  are  shifted.  And  second,  be¬ 
cause  of  the  change  in  interatomic  distance  the  u>(k)  dependence 
may  be  altered,  in  turn  changing  the  free  energy  density  field. 

In  the  presence  of  an  eccentricity  the  total  free  energy 
is  decreased  by  a  shift  of  the  Brillouin  zone  boundary  to  exclude 
a  region  of  reciprocal  space  corresponding  to  high  free  energy 
density  in  favor  of  another  region  of  lower  free  energy  density. 
This  effect  is  indicated  schematically  in  Fig.  2.  Distortion  of 
the  Brillouin  zone  boundary  corresponds  to  an  inverse  distortion 
of  the  crystal  in  real  space,  and  thence  to  magnetostriction. 

The  model  calculation  carried  out  in  Section  2  indicates  that 
the  second  effect'  (the  dependence  of  u>(if)  on  strain)  is  somewhat 
smaller  than  the  zone  boundary  effect,  and  that  both  effects  are 
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of  the  same  sign.  Both  effects  are  shown  in  Fig.  3,  where  the 
zone  boundary  effect  is  designated  as  a  "surface*  effect,  and 
the  shift  in  the  frequencies  throughout  the  zone  is  designated 
as  the  "volume'*  effect. 

As  indicated  in  Fig.  3»  the  magnetostriction  resulting  from 
a  single,  eccentricity  would  be  in  qualitative  agreement  with  the 
observations  in  iron  if  we  were  to  choose  the  eccentricity  opposite 
in  sign  and  eight  times  the  magnitude  of  that  resulting  from  dipolar 
interactions.  However,  no  plausible  source  of  that  type  of  eccen¬ 
tricity,  with  the  required  magnitude,  is  apparent  to  us.  Conse¬ 
quently,  as  an  alternative  mechanism,  we  shall  show  in  Section  3 
that  the  magnon-phonon  interaction  near  the  Brillouin  zone  bound¬ 
ary  produces  an  asymmetry  in  the  excitation  spectra  which  is  equi¬ 
valent  in  its  effect  to  a  simple  eccentricity.  The  magnon-phonon 
interaction  is  itself  sharply  temperature  dependent,  becoming 
appreciable  only  when  the  magnon-magnon  interaction  renormalizes 
the  magnon  spectrum  to  lower  the  magnon  frequencies  so  that  they 
become  degenerate  with  the  phonons  near  the  zone  edge.  This 
degeneracy  temperature  is  just  below  the  Curie  temperature  in 
iron.  The  resultant  contribution  to  \j_qq  is  therefore  very  small 
except  in  the  neighborhood  of  the  degeneracy  temperature ,' leading 
to  a  strong  maximum  just  below  the  Curie  temperature,  qualitatively 
similar  to  that  observed  in  iron. 

The  systematics  of  the  shift  and  change  of  shape  of  the  X1Qq 
maximum  with  addition  of  impurities,  and  the  criteria  for  the  pre¬ 
sence  of  the  maximum  in  other  materials,  then  follow  from  the 
relative  magnitudes  of  the  unrenormalized  magnon  and  phonon  ener¬ 
gies  in  these  materials,  and  from  the  resultant  possibility  of 
magnon-phonon  degeneracy  near  the  zone  edge. 
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II.  SIMPLE  ECCENTRICITY  MODEL 

Without  specifying  the  source  of  the  eccentricity  at  this 

stage,  and  simply  to  illustrate  the  effect  in  its  most  elementary 

form,  we  first  assume  that  the  magnon  spectrum  can  be  represented 

by  __  * 

<u(k,e)  =  u>  (k,7)  +  AYM  sin20  (l) 

o 

where  u>("k,T)  is  the  frequency  of  a  mode  of  wave  vector  £  in  an 
arbitrarily  distorted  crystal  (indicated  by  the  strain  tensor  t), 

A  is  an  undetermined  dimensionless  eccentricity  constant, 

VM  a  M  is  inserted  to  insure  that  the  eccentricity  vanish  in 

2mc 

the  symmetric  paramagnetic  state,  and  0  is  the  angle  between  ~!c  and 
ft.  This  simple  eccentricity  cannot  be  valid  at  very  low  1c,  where 
it  would  lead  to  negative  tt>(1c,7),  but  the  modes  near  the  Brillouin  , 
zone  boundary  are  of  principal  importance. 

The  dependence  of  t»o0ic,7)  on  both  7  and  7  is  taken  from  the 
Green  function  theory  of  the  Heisenberg  model.  Two  versions  of 
that  theory  have  been  given,  by  Tahir-Kheli  and  ter  Haar,^  and  by 
H.  Callen.^  Tahir-Kheli  and  ter  Haar  find  that  the  magnon  energies 
are  equal  to  simple  spin  wave  energies  multiplied  by  the  fractional 
magnetization  m(T),  The  theory  of  Callen,  while  giving  a  somewhat 
more  complicated  renormalization  which  agrees  with  the  Dyson  result 
at  low  temperature,  is  hot  radically  different  from  the  magnetization 
renormalization  at  high  temperatures.  For  simplicity,  then,  and 
because  we  are  interested  in  the  high  temperature  region,  we  simply 
adopt  the  magnetization  renormalization.  Then,  assuming  nearest 

neighbor  interaction,  the  frequency  (at  a  fixed  point  k  in  recipro- 

\  **£► 

cal  space;  becomes 


-  6  - 


hw0(k,?)  =  2JSm(T)(z  -  E  eiK’*) 


(2) 


=  2JSm(T)(z  -  E  eilt# ”°)+  2JSm(T)  E  sin  *4 


(3) 


where  6  designates  the  vectors  to  the  nearest  neighbors  (of  which 
there  are  z),  and  6q  designates  these  vectors  in  the  unstrained 
crystal.  •. 

The  free  energy  of  the  system  of  magnons  is 

F  =  E  F  =  E  £  hu>(lT)  +  p”1  E  ln(l-e-ph“^k))  (4) 

K  K  k  k 

and  we  are  interested  in  the  free  energy  difference  AF  *  F(?)- 
F(?  =  0)  in  the  strained  and  unstrained  crystal.  Furthermore  we 
are  interested  only  in  those  terms  in  AF  which  contribute  to  the 
magnetostriction  constants  *^00  an(*  ^lll»  that  they  involve 
both  the  direction  of  magnetization  (or  0)  and  the  strain^; 
denoting  these  terms  in  AF  by  AF’,  we  easily  find 

-dr  ?  (5) 


AF* 


r$  F  p  S£F 

e  B .  z . 


=h^ [$  +  n(«,T) Dftw^dT  -Ph ng  (l+n^fiWgfiw^dT  (6) 


where 


6m. 


B.Z. 


6m.  =  AyM  sin^P 

6 


■  2J  sin  R’T 


(7) 

(8) 
(9) 
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and  where  the  first  integral  is  over  the  differential  volume 
between  the  strained  and  unstrained  zones,  whereas  the  second 
integral  is  over  the  volume  of  the  unstrained  zone.  The  first 
term  in  Eq.  (6)  corresponds  to  the  surface  term  shown  in  Fig.  2. 
This  terra  has  a  zero-point  value  even  when  the  Bose  occupation 
numbers  n(u),T)  are  all  zero.  The  second  term  in  Eq.  (6)  describes 
the  effect  of  the  change  of  excitation  frequencies,  and  hence  of 
the  free  energy  density,  throughout  the  zone  when  the  crystal , is 
strained,  A  third  term,  corresponding  to  the  change  in  density 
of  states  within  the  zone,  has  been  omitted  because,  being , fully 
symmetric,  it  makes  no  contribution  to  \^qq  or  Calculation 

of  the  integrals  is  simplified  by  expanding  strains,  angular  fac¬ 
tors,  and  wave-vector  components  in  Kubic  Harmonics  and  extracting 
the  fully  symmetric  products*  We  have  carried  out  such  an  integra¬ 
tion  for  a  simple  cubic  lattice,  performing  the  numerical  integra¬ 
tions  on  the  NOL  7090*  A  check  on  the  calculations  is  provided 
by  the  exact  cancellation  of  volume  and  surface  terms  of  symmetry 
Toe*  required  by  the  invariance  of  k*^  in  the  coordinate  system 
moving  with  the  modes*  Then,  minimizing  the  sum  of  AF1  and  the 
elastic  strain  energy,  one  finds  the  equilibrium  strains,  or 
magnetostriction.  The  resultant  magnetostriction  is  shown  in 
Fig.  3,  in  which  we  show  the  portion  of  ^qq  ascribable  to  the 
surface  and  volume  terms  in  Eq,  (6).  The  temperature  dependence, 
of  m(T)  has  been  taken  from  the  empirical  magnetization  curves  of 

7 

Terry,  and  the  elastic  constants  from  the  extrapolated  data  of 

Q 

Rayne  and  Chandrasekhar.  These  are  plotted  in  Fig.  4.  It  will 

» 

be  noted  that  the  decrease  in  accounts  for  most  of 
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the  broad  \  ^  peak  on  the  simple  eccentricity  model.  The  volume 

contribution  to  is  zero  at  T  =  0  where  all  occupation'  numbers 

vanish,  and  it  falls  to  zero  again  at  T  because  of  the  renorraaliza- 
tion  of  magnon  energies  (Eq.  (8)).  The  surface  contribution  has  a 

I 

large  zero  point  contribution,  and  it  then  increases  with  increas¬ 
ing  temperature  as  the  modes  near  the  zone  boundary  become  occupied* 
The  single  adjustable  constant  appearing  in  the  results  is  the 
eccentricity  constant  A  of  Eq,  (l);  it  has  been  chosen  as  A  =  -1 6n 
in  Fig.  3,  so  that  the  theoretical  coincides  with  the  measured 

value  extrapolated  to  0°K,  For  comparison,  the  value  of  A  corre¬ 
sponding  to  the  Holstein-Primakof f^  spectrum  (i,e.,  the  dipolar- 
induced  eccentricity)  is  +2rr* 

The  corresponding  theoretical  value  of  the  shear  magneto¬ 
striction  is  properly  opposite  to  ^oq  in  sign,  but  its 

magnitude  is  only  one  third  of  the  experimental  value  at  0°K,  and 
it  drops  off  monotonically  but  too  slowly  with  increasing  temperature* 


r 
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III.  MAGNON-PHONON  INTERACTION 

We  now  consider  a  specific  mechanism  which  produces  an 

asymmetry  of  the  required  magnitude  in  the  excitation  spectrum 

of  a  ferromagnet.  Although  this  asymmetry  is  more  complex  than 

the  simple  eccentricity  considered  previously,  its  effect  is 

qualitatively  similar.  The  mechanism  to  be  considered  arises 

from  raagnon-phonon  interaction. 

The  magnon  and  phonon, spectra  each  approach  the  Brillouin 

zone  boundary  with  zero  normal  slope,  but  the  energies  vary 

considerably  over  the  face  of  the  zone.  The  phonon . frequencies 

have  been  measured  by  G.  G.  E*  Low,^  who  finds  hu>  =  0,48x10“^ 

P 

ergs  at  the  [ill]  vertex;  the  phonon  spectrum  in  this  direction 
is  shown  in  Fig*  5*  The  [ill]  vertex  is  the  principal  region  of 
low  magnon  energy  on  the  zone  surface,  and,  as  will  develop  sub¬ 
sequently,  is  therefore  of  particular  interest  to  us.  The  un¬ 
renormalized  (low  temperature)  magnon  energies,  given  by  Eq,  (2), 
can  be  obtained  sufficiently  accurately  from  the  Curie  temperature 

c  £  —13 

and  the  Green  function  theory,  ’  whence  hu>°  -  1.50x10  J  ergs  at 

m 

the  [ill]  vertex  in  iron.  The  magnon  spectrum  is  also  shown  in 
Fig.  5. 

At  low  temperatures  there  is  virtually  no  magnon-phonon 
interaction  in  iron  except  in  the  region  of  the  crossing  of  the 
two  spectra,  deep  within  the  zone.  However,  as  the  temperature 
*  increases,  the  spin  waves  energies  renormalize.  As  indicated  pre¬ 
viously,  at  high  temperature  this  renormalization  is  merely  a  t 
multiplication  of  the  spin  wave  frequencies  by  the  fractional 
magnetization  m(T).  Thus  the  magnon  spectrum  of  iron  becomes 


f 
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1 

degenerate  with  the  phonon  spectrum  at  the  zone  edge  at  that 
temperature  for  which 

m(T  )  =  <i)_/ur  (if  (10) 

p  p  m  pm 

In  actuality  the  phonon  energies  also  renormalize,  although  probably 
to  a  lesser  extent  than  the  magnons,  Accordingly,  and  in  the  ab¬ 
sence  of  a  quantitative  theory,  we  shall  neglect  this  effect. 

Because  of  the  magnetoelastic  coupling, the  phonon  and  raag- 
non  spectra  mix,  producing  new  modes  displaced  upward  and  down¬ 
ward  by  frequency  shifts  +6u>.  This  spectral  distortion  then  in- 
duces  the  zone  boundary  distortion  somewhat  as  in  the  previous 
magnon  spectrum  ellipticity  model* 

The  interaction  term  between  magnons  and  phonons  is  dictated 

by  symmetry  considerations,  and  in  the  classical  field  representa- 

10 

tion  can  be  written  in  the  form 

K  .  mC  uCj.u  d?  ■ 

mP  J  dr*  8rX  (ID 

where  the  Bfs  and  Dfs  are  phenomenological  coupling  constants* 

i 

The  effect  of  the  first  term  is  dominant  at  low  wave  vector, 
where  it  gives  rise  to  the  conventional  theory  of  magnetostriction 
and  to  the  magnon-phonon  crossing  interaction,  investigated  most 
fully  by  5chlomann.il  The  second  term  in  the  integrand  of  Eq.  (ll) 
vanishes  at  zero  wave  vector,  but  strongly  dominates  near  the 
zone  edge.  The  order  of  magnitude  of  this  term  has b een  estimated 
by  Kaganov  and  Tsukernik}^  because  it  arises  from  a  modulation  of 
the  exchange  integral  they  estimate  it  to  be  of  the  order  of  Ja^(ka)  . 
At  the  zone  boundary,  where  kaWtT,  this  would  give  a  magnon-phonon 
interaction  of  approximately  2x10“^-)  ergs/ion,  dominating  the  B  term 
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by  two  orders  of  magnitude  I  This  estimated  magnon-phonon  inter¬ 
action  would  be  as  large,  or  larger,  than  the  phonon  energy  itself, 
which  would  lead  to  an  instability  of  the  interacting  spectra,  and 
a  phase  transition*  But  the  estimate  of  Kaganov  and  Tsukernik  is 
a  very  rough  one,  and  we  shall  assume  only  that  the  magnon-phonon 
matrix  element  for  modes  near  the  zone  boundary  is  hu^,  where  u>^ 
is  less  than  but  of  the  same  order  as  the  phonon  frequency  u>p* 

The  displaced  frequencies  of  the  interacting  magnon  and 
phonon  modes  are  then  determined  by  the  secular  equation 


-  u> 


U), 


U)  -a> 
P 


=  0 


(12) 


The  two  roots  are 


s. 


•  =  < 


<0  ±  fiu> 

m  — 


(13) 


where  the  upper  sign  is  to  be  taken  if  wra>*p  and  the  lower  sign 

if  and  where 

ra  p 9 


2,i 


-  -  *l<V<“pl  +  +4<"i  ] 


(14) 


Thus  the  displacement  of  the  modes  reaches  a  maximum  value  of 
when  u>_  =  <u ,  and  it  falls  to  ^  when  u>  -u»  =  + 

The  interaction  matrix  element  is  certainly  strongly 
dependent  on  the  angle  between  the  wave  vector  of  the  modes  and 
the  magnetization,  as  has  been  shown  at  the  cross-over  region  at 
low  1c  by  Schlomann.^  It  is  obvious  by  symmetry  that  a  spin  wave 
propagating  along  the  magnetization  cannot  couple  to  a  longitudinal 
phonon.  However,  the  coupling  can  be  large  for  modes  propagating 
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perpendicular  to  the  magnetization  direction.  We  thus  summarize 
•the  above  considerations  by  assuming  that  0<i>(0)  is  a  function  of 
the  angle  0  between  k  and  M; 

6U) ( 0 )  =  g(»ra-wp)  Sin  0  (15) 

where  g  is  of  the  order  of  dip,  falling  to  half  its  value  when 

u>  -id  ~  +  .2.e  ( 0 ) . 
m  p  —  28 

Again  the  free  energy  of  the  crystal  is  decreased  by  a 
distortion  of  the  Brillouin  zone,  excluding  regions  of  small 
splitting  in  favor  of  regions  of  large  splitting.  As  in  Eq.  (6) 
the  surface  contribution  to  the  change  in  free  energy  is 


AF" 


=  /  [F  (<#  +4w)+F  (u>  +6tu)]dT 
<J  k  m  X  p 


(16) 


At  the  degeneracy  temperature  this  change  in  free  energy  becomes 


AF"  =  -ph2  /  n*  (l+nj-)  g2  sin20  dT 

■  ,J*  K 


(17) 


Comparison  of  this  equation  with  the  surface  term  in  Eq.  (6)  shows 

that  the  same  value  of  magnetostriction  at  will  be  obtained  if 
2 

g  and  A  are  related  by 


A  =  -»p  _ 

[l-exp(-p^h»p)]vM0d)p  (18) 

Consequently  our  previous  result  that  a  value  of  A  =  -16tt  was 
required  to  fit  the  value  of  ^qq  at  implies  that  this  fit 
will  be  obtained  with  hg  -  10“^  ergs.  This  required  value 
actually  is  considerably  less  than  the  value  estimated  from  the 
magnon-phonon  matrix  element.  1 

We  first  estimate  the  position  and  shape  bf  the  peak  in 
pure  irr-n,  and  then  discuss  the  effect  of  silicon  alloying. 
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The  peak  temperature,  from  Eq.  (10),  and  Fig.  5,  is  that  tem¬ 
perature  at  which  m(T  )  =  From  the  experimental  raagnetiza- 

tion  curve  of  Fig.  4  we  find  Tp  =  0.96  Tc  =  1000°K.  This  figure, 
somewhat  higher  than  the  experimental  value,  is  reduced  by  the 
observation  that  both  the  phonon  and  magnon  frequencies  vary 
somewhat  over  the  zone  surface.  Hence  magnon-phonon  degeneracy 
is  reached  at  different  temperatures  at  different  points  on  the 
surface*  The  observed  maximum  therefore  should  be  broadened  and 
depressed  from  our  initial  theoretical  estimate. 

In  addition  to  this  " geometrical11  broadening,  at  least 
two  other  effects  broaden  the  magnetostriction  maximum.  At  high 
temperatures  the  lifetimes  of  the  magnons  and  phonons  become 
short,  so  that  the  excitations  are  themselves  broadened  consider¬ 
ably.  Furthermore  there  is  an  intrinsic  breadth  of  the  maximum 
which  arises  from  the  fact  that  the  magnons  and  phonons  interact 
even  at  temperatures  other  than  T^. 

We  recall  that  6u>  is  a  maximum  at  T  ,  and  that  it  falls  to 

P 

half  its  maximum  at  the  temperatures  T  i  determined  by 


±i 

"(T±i> '  -  "‘v*  i  £  • 

m  m 


(18) 


If  we  assume  arbitrarily  that  m^-Wp,  then  m(T_^)  =*  0.8  and 
m(T+^)  s  0,  whence  T_^  is  roughly  800°K  and  of  course 
T^  =  Tq  =  1043°K.  Like  the  peak  temperature,  these  "half  width" 
points  are  again  too  high,  but  the  theory  predicts  correctly  that 
the  line  will  be  broadened  asymmetrically  because  of  the  shape  of 
the  magnetization  curve,  with  the  low  temperature  rise  flitter 
than  the  sharp  high  temperature  drop-off. 
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The  magnon-phonon  interaction  mechanism  also  accounts 
naturally  for  the  considerable  effect  of  small  additions  of 
silicon.  The  position  of  the  magneto-striction  peak  is  de¬ 
termined  by  the  phonon  modes  at  the  zone  boundary,  and  these 
modes  are  radically  perturbed  by  the  addition  of  light  im¬ 
purities.  In  fact  the  tendency  of  such  impurities  is  to 
split  localized  modes  off  the  top  of  the  phonon  band;  even 

if  such  localization  is  not  accomplished,  the  modes  of  short 
m 

wave-length  (\  ot  a)  are  shifted  upward  in  frequency.  Conse¬ 
quently,  the  addition  of  silicon  shifts  the  magnetostriction 
peak  to  lower  temperatures,  and  broadens  it  by  moving  the 
peak  to  a  region  of  smaller  ^  J  •  These  effects  appear  in 

Fig.  1. 


If  is  less  than  «b^,  rather  than  greater  as  in  iron, 

no  degeneracy  temperature  occurs.  Magnon-phonon  interaction 

will  then  give  a  contribution  to  the  magnetostriction  at  zero 

temperature  if  <#p-u>°,£(i>^ ,  and  this  contribution  will  decrease 

with  increasing  temperature.  This  appears  to  be  the  case  in 

nickel.  Hard  materials  (high  0)^)  with  low  Curie  temperatures 

should  not  show  a  high-temperature  magnetostriction  maximum. 

The  optimum  condition  for  the  observation  of  the  high- 

temperature  magnetostriction  maximum  is  that  be  perhaps  twice 

to  five  times  as  large  as  <#  .  For  if  tu°  is  enormously  larger 

p  m 

than  #>p  the  degeneracy  temperature  becomes  practically  identical 
to  the  Curie  temperature,  and  the  peak  becomes  very  narrow.  The 
situation  in  iron  is  apparently  very  close  to  the  optimum 
conditions. 
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Fig.  la. 

Fig.  lb. 

Fig.  2. 


Fig.  3. 


Fig.  4. 


FIGURE  CAPTIONS 

Magnetostriction  constant  ^qq  vs  temperature,  of  iron 
and  silicon  iron,  according  to  E,  Tatsumoto  and  T.  Okaraoto. 
J.  Phys.  Soc.  Japan  1588  (1959). 

Magnetistriction  constant  A  vs  temperature,  of  iron 
and  silicon  iron,  according  to  Tatsumoto  and  Okamoto. 

Zone  boundary  distortion  effect. 

An  ellipticity  in  the  magnon  spectrum  with  respect  to  the 
magnetization  direction  creates  an  asymmetry  in  the  free 
energy  density  field.  The  free  energy  is  then  lowered  by 
extending  the  zone  boundaries  to  include  regions  of  low 
free  energy  and  density  and  contracting  the  boundaries  to 
exclude  the  high  energy  regions. 

Magnetostriction  constant  ^qq  as  a  function  of  tempera¬ 
ture  arising  from  the  simple  eccentricity  model.  The 
ellipticity  coefficient  has  been  adjusted  to  reproduce 
the  experimental  magnetostriction  extrapolated  to  0°K. 

The  "surface"  contribution,  as  indicated  in  Fig.  2,  and 
the  "volume"  term,  due  to  the  change  in  energy  of  the 
magnons  throughout  the  zone  when  the  crystal  is  strained, 
are. shown  separately,  as  is  their  sum  and  the  experimental 
curve  of  Fig.  1. 

The  experimental  reduced  magnetization  of  iron  as  a  func¬ 
tion  of  temperature,  according  to  E.  M.  Terry,  Phys.  Rev. 

133  (1910).  The  elastic  stiffness  constants  c^(T)  ”  Cj^T) 
and  c^(T)  of  iron  have  been  measured  by  J.  A.  Rayne  and 
B.  S.  Chandrasekhar,  Phys.  Rev.  122.  1714  (1961)  frdm 
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FIGURE  CAPTIONS  -  Cont. 

4.2°K  to  500°K.  K.  Nakamura,  Sci.  Rept.  Tohoku  Univ.  2*>. 
3 64  (1955-57)  measured  the  compliance  constants  of  iron 
in  the  range  295°K<T<795°K.  The  stiffness  constants  em¬ 
ployed  by  us  are  the  data  of  Rayne  and  Chandrasekhar 
extrapolated  parallel  to  the  measurements  of  Nakamura 
( converted ) . 

Fig,  5.  Phonon  energy  is  a  function  of  propagation  vector  in  the 
[ill]  direction  in  iron,  according  to  G,  G,  E,  Low,  Proc. 
Phys.  Soc.  22,»  479  (1962).  The  unrenormalized  (0°K) 
magnon  energy  vs  propagation  vector  in  the  same  direction 
in  reciprocal  space  is  also  shown. 
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MAGNETOSTRICTION  CONSTANT  x 


REDUCED  MAGNETIZATION  m(T) 


ELASTIC  CONSTANT (DYNES/cm2) 5 


